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ABSTRACT 


This  report  describes  the  preliminary  evaluation  of  the  Iranian 
Long  Period  Array  (ILPA).  This  evaluation  was  performed  by  Texas  Instru- 
ments Incorporated  at  the  Seismic  Data  Analysis  Center  in  Alexandria, 
Virginia. 

The  major  areas  of  study  in  this  evaluation  are: 

• Evaluation  of  the  data  quality  and  sources  of  data  errors 

• Estimation  of  beamforming  gains 

• Estimation  of  seismic  event  detection  thresholds 

• Determination  of  seismic  event  Mg-m^  relationships. 

Conclusions  regarding  the  above  points  and  plans  for  future 
work  are  also  presented. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and  this 
document  is  subject  to  later  revision  as  may  be  necessary.  The  views  and  con- 
clusions presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 
INTRODUCTION 


This  report  presents  the  results  of  a preliminary  evaluation  of 
the  seven  element  Iranian  Long-Period  Array  (ILPA).  Since  data  did  not  begin 
to  arrive  until  the  end  of  May  1976,  it  was  necessary  to  sharply  curtail  the 
goals  of  this  evaluation.  In  the  limited  time  remaining  in  the  contract  period, 
emphasis  was  placed  on  evaluating  the  quality  of  the  ILPA  data  and  obtaining 
estimates  of  ILPA  detection  capability.  The  specific  areas  of  investigation 
include: 

• Data  quality 

• Sources  of  data  errors 


e Beamforming  gains 

• Seismic  event  detection  thresholds 


• Seismic  event  M^m^  relationships. 

A brief  description  of  the  Iranian  Long  Period  Array  iB  given 
in  Section  II.  The  data  base  and  data  processing  methods  are  described  in 
Section  III.  Data  quality  and  sources  of  data  errors  are  discussed  in  Section 
IV.  Beamforming  gains,  detection  thresholds,  and  M -m,  relationships  are 

8 D 

discussed  in  Section  V.  Finally,  the  results  of  this  preliminary  evaluation  of 
the  ILPA  data  are  summarized  in  Section  VI. 
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SECTION  II 

IRANIAN  LONG-PERIOD  ARRAY  SYSTEM  DESIGN 


A.  DESCRIPTION  OF  THE  ARRAY 

The  Iranian  Long  Period  Array  is  a seismic  installation  com- 
prised of  a central  recording  station  and  an  array  of  seven  remote  sites.  The 
locations  of  these  sites  are  listed  in  Table  II- 1 and  shown  in  Figure  II- 1. 

Each  remote  site  is  composed  of  a three -component  sensor 
subsystem,  a data  acquisition  subsystem,  a telemetry  subsystem,  and  a 
power  subsystem.  The  three-component  seismometer  of  the  sensor  subsys- 
tem is  located  in  a 100  meter  borehole  to  reduce  wind-generated  noise.  The 
seismometer,  a signal  conditioning  unit,  and  a remote  control  unit  make  up 
the  sensor  subsystem.  The  data  acquisition  subsystem,  housed  in  a small 
building  near  the  borehole,  contains  a six-channel  multiplexer,  a sample - 
and-hold  module,  a binary  gain-ranging  amplifier,  an  analog-to-digital  con- 
verter, and  data  formatting  and  control  logic.  The  data  acquisition  subsystem 
converts  signals  from  the  sensor  subsystem  to  the  correct  format  for  trans- 
mission to  the  central  recording  station. 

The  ri  mote  site  telemetry  subsystem  transmits  signals  from 
the  data  acquisition  subsystem  to  the  central  recording  station  and  receives 
command  signals  from  the  central  recording  station.  The  equipment  at  each 
site  consists  of  an  antenna,  a duplexer  which  permits  transmitting  and  receiv- 
ing with  a single  antenna,  a transmitter,  a receiver,  transmit  and  receive 
modules,  and  remote  control  logic.  One  telemetry  subsystem  includes  a 
relay  station,  located  at  site  5,  which  relays  data  between  site  6 and  the 
central  recording  station.  This  relay  is  necessary  because  site  6 does  not 
have  line-of- sight  with  the  central  recording  station. 
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TABLE  II- 1 j 

1 I 

REMOTE  SITE  COORDINATES  ; 


Site 

Location 

Distance  From  Reference  Site 
(km) 

1151 

H 

North 

East 

1 (ref) 

35°24'58.  3" 

50°41 1 1 9.  5" 

0.  0 

0.  0 

2 

35°39'46.  1" 

50°53'51. 5" 

27. 277 

19.  035 

3 

35°28'34.  0" 

51°01 !25.  5" 

6.  217 

30.377 

4 

35°14'19.  3" 

50°54'04. 1" 

-19.536 

19.  162 

5 

35°1 2'46.  2" 

50°34'52.  0" 

-22.415 

-9.  830 

6 

35°28'25.  2" 

50°25'32.  2" 

5.815 

-23. 775 

7 

35°42' 10.  1" 

50°36'32.  0" 

31.  700 

-6.  951 

The  central  recording  station  processes  and  records  data  re- 
ceived from  each  site  of  the  array.  It  contains  the  equipment  necessary  to 
transmit  commands,  receive  and  process  data  from  the  remote  sites,  and 
record  the  results  of  this  processing  on  visual  and  magnetic  recorders.  A 
timing  subsystem  for  synchronization  of  operations  at  the  seven  remote  sites 
and  the  central  recording  station  is  also  part  of  the  central  recording  station. 
The  equipment  housed  in  the  central  recording  station  are  the  station  processor, 
the  visual  recording  subsystem,  the  magnetic  digital  tape  recording  subsystem, 
the  timing  subsystem,  the  telemetry  subsystem,  and  the  power  subsystem. 

The  station  processor  contains  the  computing  and  control  ele- 
ments of  the  system,  supplied  by  the  central  processing  units.  Each  of  the 
central  processing  units  contains  the  same  program  and  is  capable  of  perform- 
ing all  on-line  and  off-line  requirements  of  the  system.  The  central  process- 
ing unit  which  is  manually  placed  on-line  assumes  the  responsibility  of  per- 
forming on-line  operations  and  transmits  to  the  off-line  central  processing 
unit  at  one -second  intervals.  This  transmission  suppresses  on-line  opera- 
tions by  the  off-line  central  processing  unit.  Failure  to  receive  a transmis- 
sion activates  the  on-line  coding  of  the  off-line  central  processing  unit, 
switching  the  roles  of  the  central  processing  units. 

A teleprinter  interfaces  with  each  of  the  central  processing 
units.  These  provide  man -machine  communications.  Commands  are  given 
to  the  system  by  keyboard  entry  and  the  system  status  is  printed  out  automati- 
cally. 

The  visual  recording  system  converts  the  digital  signals  from 
the  remote  sites  back  to  analog  form  for  display  on  drum  recorders  and  devel- 
ocorders.  The  develocorder a record  the  data  on  film,  process  the  data,  and 
store  it.  Viewing  screens  on  the  develocorder s permit  viewing  of  the  film. 

A second  visual  display  is  made  by  the  drum  recorders  on  30x90  centimete.. 
recording  paper. 
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The  digital  magnetic  recording  system  records  three-comnonent 
long-period  data  from  the  remote  sites.  One  recorder  is  reserved  as  ba>.  ,^up 
in  case  of  an  on-line  recording  drive  failure.  This  system  is  also  used  for 
the  tasks  of  providing  data  for  beamforming  and  display  and  editing  data  to 
another  tape,  It  was  intended  that  data  would  be  transmitted  to  the  United 
States  by  satellite  relay.  To  accomplish  this,  the  data  is  reformatted  into  the 
satellite  format  of  Figure  II-2.  To  date,  the  data  is  not  relayed  by  satellite 
but  is  recorded  on  magnetic  tape  in  this  format  at  1LPA.  The  tapes  are  then 
shipped  to  the  United  States. 

The  telemetry  system  transmits  command  and  control  data  to 
the  remote  sites  as  well  as  calibration  signals  to  the  sensor  Bystems.  This 
system  also  receives  data  from  the  remote  sites  via  seven  receivers.  Each 
receiver  drives  a receiver  module  which  interfaces  through  a separate  I/O 
port  to  the  central  processing  units.  These  interfaces  perform  Berial-to- 
parallel  conversion  of  incoming  data  and  present  each  central  processing  unit 
with  data  on  an  interrupt  basis. 

The  power  system  is  designed  to  provide  uninterrupted  con- 
trolled power  to  the  central  processing  station.  When  this  system  is  com- 
pleted, power  from  a motor  generator  will  be  automatically  substituted  via  a 
transfer  switch  for  commercial  power  during  commercial  power  failures. 

The  generator  will  be  capable  of  sustaining  power  for  eight  hours  without  re- 
fueling. The  transfer  switch  will  allow  for  a delay  of  up  to  five  minutes  be- 
fore the  generator  is.  started.  During  this  delay,  battery  power  will  maintain 
the  system.  These  batteries  will  power  the  system  for  one  hour. 

B.  THE  DATA  RECORDING  FORMAT 

The  primary  output  of  the  data  recording  system  available  to 
us  is  the  digital  magnetic  tape.  This  tape  is  nine-track,  recorded  In  the  satel- 
lite tape  format  at  800  bits  per  inch,  using  two's  complement  binary  arithmetic. 
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Satellite  Tape  Format 


Format  of  Satellite  Output  Buffer* 


Bytes 

0-1 

F09A16 

Bytes 

2-3 

'IL' 

Bytes 

4-11 

Status  Bytes 

Bytes 

12-17 

Time  Code 

Bytes 

18-59 

Long -Period 
Data 

Bytes 

60-99 

Short-Period 

Data 

Bytes 

100-101 

C8C816 

Each  record  of  20  buffera 
contains  20  seconds  of  data. 


24  hours  per  tape 
800  bpi 
9 track  tape 
8.  5,!  reels 


* Literals  are  in  EBCDIC* 


FIGURE  II- 2 


Sync 

Station 

Code 


End  of 
Frame 
Indicator 


SATELLITE  TAPE  FORMAT 
(PAGE  1 OF  3) 
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Format  of  Status  Bytes 


Byte  0 

Status  Bits 
Site  1 

Byte  1 

Status  Bits 
Site  2 

Byte  2 

Status  Bits 
Site  3 

Byte  3 

Status  Bits 
Site  4 

Byte  4 

Status  Bits 
Site  5 

Byte  5 

Status  Bits 
Site  6 

Byte  6 

Status  Bits 
Site  7 

Byte  7 

Reserved 

Format  of  the  Time  Code* 

Bytes  0-1 

OYYD 

Bytes  2-3 

DDHH 

Bytes  4-5 

• MMSS 

Bit 

0 

Sync  Error 

Bit 

1 

Calibration  in 
progress 

Bit 

2 

Deleted  from 
beamforming 

Bit 

3 

Faulty  or  Missing 
Long-Period  Data 

Bit 

4 

Faulty  or  Missing 
Short-Period  Data 

Bit 

5-7 

Reserved 

(The  above  bits  are  set  to 
one  if  the  condition  existB 
or  to  zero  if  the  condition 
does  not  exist. ) 


Year:  0<YY<9910 

Julian  Day:  1<DDDS36510 

(on  leap  year, 
1<  DDDS366i0) 

Hour:  0<  HH<L23j0 

Minute:  0<MM^591g 

Second:  0£SSS59^q 


>'.<  Contents  are  in  packed  BCD. 
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SATELLITE  TAPE  FORMAT 
(PAGE  2 OF  3) 
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Format  of  Long-Period  Data 


Bytes  0-1 

Site  1,  Vertical  Component 

Bytes  2-3 

Site  1,  N-S  Component 

Bytes  4-5 

Site  1,  E-W  Component 

Bytes  6-7 

Site  2,  Vertical  Component 

Bytes  8-9 

Site  2,  N-S  Component 

Bytes  10-11 

Site  2,  E-W  Component 

( 

• 

Bytes  36-37 

Site  7,  Vertical  Component 

Bytes  38-39 

Site  7,  N-S  Component 

Bytes  40-41 

Site  7,  E-W  Component 

Each  component  is  in 
gain -ranged  format. 


Format  of  Short-Period  Data 


Bytes  0-1 
Bytes  2-3 


Bytes  38-39 


Channel  Datum  for  Zeroth 
Short-Period  Frame 


Channel  Datum  for  First 
Short-Period  Frame 


Channel  Datum  for  Nine- 
teenth Short-Period  Frame 


Each  datum  is  in 
gain -ranged  format. 


(Only  one  channel  of 
short-period  data  can 
be  placed  in  the  above 
buffer  at  any  one  time. ) 


FIGURE  II -2 

SATELLITE  TAPE  FORMAT 
(PAGE  3 OF  3) 
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The  satellite  tape  format  is  shown  in  Figure  II-2,  which  illustrates  the  encod- 
ing of  seven  sites  of  three -component  long-period  data  and  one  site  of  vertical- 
component.  short-period  data. 

The  ILPA  data  word  format  consists  of  two  bytes  containing 
eight  bits  each  in  the  following  format: 

• Byte  1:  The  four  most  significant  bits  define  the  gain  factor. 

The  next  bit  is  the  sign  bit.  The  last  three  bits  are  the 
three  most  significant  bits  of  the  mantissa. 

• Byte  2:  These  eight  bits  are  the  eight  least  significant  bits  of 

the  mantissa.  The  mantissa  is  encoded  in  two's  com- 
plement binary  arithmetic. 

The  voltage  registered  by  a sensor  may  be  calculated  by  the 

formula: 


Voltage 


Mantissa  * 2 


“(gain  code) 


7FF 


* 10 


16 


10  volts' 


The  gain  code  varies  fromO  to  4 for  short-period  data  and 
from  0 to  for  long-period  data.  The  number  10^  is  the  full-scale  input 
signal  accepted  by  the  analog-to-digital  converter.  Using  the  above  formula, 
the  smallest  signal  which  may  be  resolved  is: 


Long-Period:'  VMIN 


1 » 2~A16 

7FFlT" 


* 10 


10 


4.  773  microvolts  ( |jlv) 


Short-Period:  V. 

MIN 


1 » 


7FF 


* 10 


16 


10 


0.  30  5 millivolts  (mv) 


Thus,  for  long-period  data,  one  computer  count  equals  4.  773|uv,  while  for 
short-period  data,  one  computer  count  equals  0,  305  mv. 
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The  sensor  filters  are  set  for  output  as  follows; 


volts 

Long-Period;  0*  1 ^ , ■ at  25  seconds  period 


micron 


Short-Period;  5.  0 -v°^ts — at  i second  period  . 

micron  r 

Therefore,  the  smallest  signal  in  millimicrons  which  may  be  resolved  is: 


Long-Period:  SIG 


4.  773  >1*  10  volts  * 103  milUmicrons/micron 


MIN  0.  1 volts/micron 


= 0.0477  3 millimicrons  (m\i) 


Short-Period:  SIG 


MIN 


. -3 


0.  30  5 * 10  volts  * ,„3  . 

■s— r-  — r-— i — ; * 10  millimicrons/ micron 

5.0  volts/micron 


= 0.061  millimicrons  (m(j.) 


From  this  we  see  that  the  data  is  quantized  at  one  computer  count  per  0.04773 
millimicrons  (20.951  cc/ xa.fi)  for  long-period  data  at  25  seconds  period.  The 
short-period  data  is  quantized  at  one  computer  count  per  0.061  millimicrons 
(16.  393  cc/m/i)  at  one  second  period. 

A more  detailed  description  of  the  ILPA  system  can  be  found 
in  the  operation  and  maintenance  manual  for  the  ILPA  Beismic  system. 
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SECTION  III 

DATA  BASE  AND  METHOD  OF  ANALYSIS 


A.  DATA  BASE 

For  our  first  look  at  the  Iranian  Long-Period  Array  data,  we 
selected  the  data  recorded  during  the  month  of  May  1976  since  this  was  the 
earliest  data  we  received.  All  events  listed  in  the  Norwegian  Seismic  Array 
(NORSAR)  event  lists  for  May  1976  which  had  epicenters  in  or  near  the  Euras- 
ian landmass  were  selected  for  processing.  No  effort  was  made  to  eliminate 
events  which,  based  on  origin  time  and  epicenter  location,  could  be  expected 
to  arrive  at  ILPA  in  the  same  time  gate  ("mixed  signals").  This  gave  us  a 
total  of  281  events  with  epicentral  distances  ranging  from  0.6  to  74.0  degrees. 
Figure  III- 1 shows  that  this  event  population  can  be  broken  down  into  a near- 
field  subset  (0°  to  20°),  a near  teleseismic  subset  (20°  to  50°), and  a far  tele- 
seismic  subset. 

We  note  that  at  the  time  the  data  base  was  formed,  we  had  no 

depth  information  for  these  events.  Thus,  it  is  entirely  possible  that  some  of 

these  events  occurred  at  depths  significantly  greater  than  33  km  (normal  depth). 

Our  interest  in  the  depth  of  occurrence  lies  in  the  fact  that  deep  events  have 

significantly  lower  surface -wave  magnitudes  (M^)  than  do  shallow  events. 

Thus,  a set  of  deep  and  shallow  events  will  have  a greater  M variance  than 

s 

will  a set  of  shallow  events  akne.  No  presumed  explosions  were  known  by  us 
to  have  occurred  during  May;  so  we  cannot  discuss  the  question  of  earthquake- 
presumed  explosion  discrimination.  The  seismic  parameters  of  the  data  base 
are  listed  in  Appendix  A. 
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The  noise  data  base  was  formed  by  searching  the  NORSAR  event 
lists  for  time  intervals  of  at  least  one  hour  in  duration  for  each  day  during 
which  no  signals  would  be  expected  to  arrive  at  ILPA.  Segments  4096  seconds 
in  length  were  processed  and  visually  examined  for  unreported  signals  or  sys- 
tem malfunctions.  If  either  was  observed,  a new  time  interval  for  that  day 
was  selected  . 

B.  METHOD  OF  DATA  PROCESSING  AND  ANALYSIS 

Five  basic  programs  were  used  to  perform  the  data  processing. 
The  first  of  these  carried  out  a check  for  incorrect  record  lengths  and  a prob- 
lem we  term  the  "data -shift  error,  " where  the  data  is  shifted  one  byte  out  of 
the  correct  format  positions.  This  gave  us  an  idea  of  the  quality  of  the  data 
on  each  tape.  With  this  information  we  could  avoid  wasting  time  trying  to 
process  data  from  unreadable  tapes. 

The  second  program  was  also  used  to  chock  the  data  quality. 
This  program  checked  for  parity  errors  and  timing  errors,  printing  out  the 
type  of  error  and  the  time  at  which  it  occurred. 

The  third  program  edited  the  desired  time  segments  from  the 
data  tapes,  performed  quality  checks,  computed  the  trace  means  for  the 
three  components  of  motion  of  each  site,  and  output  these  time  segments  to  a 
permanent  hold  tape.  The  quality  check  part  of  this  program  printed  out  mes- 
sages indicating  the  presence  of  parity  errors,  timing  errors,  clipped  data, 
and  spiked  data.  Also  printed  out  is  a summary  of  segment  powers,  which 
can  be  UBed  to  find  bad  sites.  (Bad  sites  are  those  sites  which  are  dead, 
contain  uncorrectable  spikes,  or  display  abnormally  high  power  levels.  ) 

This  edit  program  was  created  b}  adapting  the  search  and 
read  subroutines  of  the  Alaskan  Long  Period  Array  (ALPA)  edit  program  to 
handle  the  ILPA  satellite  tape  format.  The  ALPA  array  parameters  in  the 
program  were  replaced  with  those  for  ILPA.  By  adapting  this  program  to 
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handle  ILPA  data,  all  succeeding  programs  of  the  array  processing  package 
could  be  applied  to  the  ILPA  data  edits  without  modification.  The  modified 
edit  program  was  checked  by  manually  decoding  hexadecimal  dumps  of  por- 
tions of  an  ILPA  tape  and  comparing  the  results  with  the  output  of  the  edit 
program. 

The  fourth  program  performed  trace  mean  removal,  rotation 
of  the  three  components  of  motion  from  their  recorded  vertical,  north,  east 
(V,  N,E)  configuration  to  a vertical,  transverse,  radial  (V,  T,  R)  configuration, 
and  beamforming  of  the  data.  By  rotating  the  data,  we  separate  the  surface 
waves  recorded  on  the  horizontal  components,  resulting  in  two  components  of 
Rayleigh  motion  (V  and  R)  and  one  component  of  Love  motion  (T),  Noise 
samples  retained  their  V,  N,  E configuration. 

The  beamforming  process  performed  by  this  program  increases 
the  signal-to-noise  ratio  by  suppressing  random  noise  by  a factor  approxi- 
mately equal  to  the  square  root  of  the  number  of  sites  used.  Three  compon- 
ents of  motion  from  a reference  site  (normally,  site  1)  and  the  three  beams 
are  output  to  a permanent  hold  tape  with  appropriate  annotation.  This  pro- 
gram was  checked  by  manually  decoding  a hexadecimal  dump  of  a portion  of  a 
data  tape  and  hand -computing  the  beam  output,  This  was  then  compared  to 
the  output  of  the  beamforming  program. 

The  fifth  program  bandpass  filtered  the  reference  site  and  beam 
traces  and  output  plots  suitable  for  analysis. 

The  data  processing  method  is  illustrated  in  Figure  III-Z. 

Once  the  data  were  plotted,  signal  analysis  consisted  of  deter- 
mining whether  the  event  was  detected  on  the  single-site  and  beam  traces  and 
measuring  amplitude  and  period  information  on  detected  events  for  surface- 
wave  magnitude  computations.  Signal-to-noise  ratios  for  both  single-site  and 
beam  data  were  measured  on  the  bandpass-filtered  traces  to  allow  estimation 
of  beamforming  gains.  All  plots  were  visually  checked  for  malfunctions  and 
mixed  signals. 
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SECTION  IV 
DATA  QUALITY 


A.  DATA  TAPE  ERRORS 

When  the  ILPA  edit  program  was  first  created,  it  was  set  up 
to  read  data  in  the  satellite  format  exactly  as  shown  in  Figure  II-2  of  Section 
II.  Two  problems  appeared  immediately  when  the  program  was  tested  on 
ILPA  data.  First,  it  was  found  that  each  tape  has  two  file  marks,  at  the  start 
of  the  tape  which  are  not  mentioned  in  the  format  description.  Once  we  were 
aware  of  this  problem,  it  was  easy  to  correct  by  adding  a statement  in  the 
program  to  skip  over  these  file  marks  before  beginning  to  read  the  data. 

The  second  problem,  which  we  may  term  the  "data-shift" 
error,  is  illustrated  by  Figure  IV-1.  (In  previous  discussions  of  ILPA  data, 
we  have  referred  to  this  as  the  "sync  error.  " We  here  change  the  name  to 
avoid  confusion  with  a transmission  problem  at  ILPA  which  is  also  called  the 
"sync  error.  ")  When  this  data-shift  error  occurs,  we  find  that  all  the  data 
has  been  shifted  down  one  bvte  in  the  format.  Each  data  record,  according  to 
the  satellite  data  format,  is  to  start  with  the  hexadecimal  characters  FO  and 
9A,  followed  by  C9,  D3  ('IL1)  and  ending  with  the  characters  C8,  C8.  When 
the  data-shift  error  occurs,  the  first  hexadecimal  characters  are  C8,  FO. 

The  9A  character  is  in  the  location  reserved  for  the  C9  character  (I  of  'IL') 
and  so  on  through  the  record.  This  problem  occurs  frequently.  Of  the  31 
tapes  containing  May  1976  data,  two  tapes  could  not  be  read  and  the  remaining 
29  showed  this  data-shift  error  for  more  than  half  of  the  time  period  recorded. 
It  was  found  that  once  this  shift  occurred,  it  continued  for  many  hours. 
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THE  DATA -SHIFT  ERRO 


The  source  of  this  problem  is  not  clear.  'We  sent  a copy  of 
one  of  the  ILPA  satellite  format  data  tapes  to  Mr.  Wayne  Ellis  of  Texas  In- 
struments Incorporated  in  Dallas  with  a list  of  the  times  at  which  we  found 
these  data-shift  errors.  Using  hexadecimal  dumps  of  this  tape,  he  was  unable 
to  find  these  errors.  Furthermore,  during  a subsequent  trip  to  ILPA,  he  was 
unable  to  find  a cause  for  these  errors.  After  discussing  this  problem  with  an 
IBM  computer  engineer  in  Dallas,  he  suggested  that  the  source  of  these  errors 
may  lie  in  the  tape  transport  of  the  IBM  360/44  located  here,  causing  the  com- 
puter to  misread  the  inter-record  gaps. 

Whatever  the  cause  of  this  problem,  we  still  had  to  process  the 
data.  Therefore,  the  ILPA  edit  program  was  altered  so  that  it  firBt  searched 
for  the  hexadecimal  characters  C9,D3  ('IL')  and  aligned  the  data  relative  to 
their  location.  Thus,  the  edit  program  now  sees  no  difference  between  cor- 
rectly formatted  data  and  data  containing  the  data-shift  error. 

Once  we  began  processing  the  events  of  our  data  base,  we  began 
encountering  parity  errors  and  timing  errors.  Using  the  second  data  process- 
ing program  described  in  Section  III,  we  counted  the  parity  and  timing  errors 
on  each  readable  tape.  The  results  are  listed  in  Table  IV-1.  (Some  of  the 
tapes  could  not  be  read  by  this  program  due  to  excessive  parity  or  timing 
errors  or  missing  end-of-tape  marks.)  Upon  encountering  a parity  error, 
the  edit  program  prints  out  an  error  message  indicating  the  presence  and 
approximate  location  of  the  error.  It  then  backspaces  the  record  and  attempts 
to  read  the  data  a second  time.  The  data  is  accepted  even  if  the  parity  error 
occurs  on  the  second  read.  The  effect  of  these  parity  errors  on  the  data  does 
not  appear  to  be  severe.  When  comparing  plots  of  the  data  with  the  parity 
error  listing  generated  by  the  edit  program,  with  rare  exceptions  we  do  not 
find  that  the  parity  error  had  any  effect  oh  the  data.  When  the  parity  error 
occurred  in  conjunction  with  a timing  error,  the  plot  showed  a spike  in  the 
data. 
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TABLE  IV- 1 

SUMMARY  OF  PARITY  AND  TIMING  ERRORS 


Tape 

Julian  Date 

Number  of 
Parity  Errors 

Number  of 
Timing  Errors 

1 

122-123 

10 

4 

2 

122-124 

4 

4 

3 

124-125 

B 

8 

4 

125-126 

27 

8 

5 

126-127 

0 

11 

6 

127-128 

157 

7 

8 

129-130 

1 

2 

10 

131-132 

4 

7 

11 

132-133 

4 

13 

12 

133-134 

0 

8 

14 

135-136 

5 

2 

17 

138-139 

2 

13 

18 

139-140 

0 

8 

19 

140-141 

8 

10 

20 

141-142 

2 

18 

21 

142-143 

4 

13 

22 

143-144 

0 

8 

2 3 

144-145 

1 

18 

24 

145-146 

2 

9 

25 

146-147 

0 

5 

26 

147-148 

135 

2 

27 

148-149 

0 

3 

28 

149-150 

0 

5 
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Upon  encountering  a timing  error  (time  of  sample  i + 1 is  not 
equal  to  time  of  sample  i plus  the  sample  rate),  the  edit  program  prints  out 
a message  indicating  the  presence  and  location  of  the  error,  accepts  the  data 
associated  with  the  timing  error,  and  continues.  As  long  as  there  is  no  tim- 
ing error  at  the  time  the  edit  is  to  start,  the  edit  will  run,  since  timing  is 
maintained  in  the  program  by  a reference  timing  word.  Timing  errors  for 
data  points  following  the  first  point  are  detected  by  comparing  the  time  on 
tape  associated  with  each  data  point  with  the  reference  timing  word.  We  gen- 
erally see  a spike  on  plots  of  the  data  when  a timing  word  error  occurs.  This 
implies  that  the  entire  data  sample  is  in  error  and  not  just  the  timing  word. 

B.  NUMBER  OF  GOOD  SITES 

One  indication  of  the  quality  of  the  ILPA  data  1b  the  number  of 
sites  considered  to  be  acceptable  for  beamforming.  A site  is  rejected  from 
beamforming  if  any  of  the  following  occurs: 

• The  site  is  dead. 

• One  or  more  components  of  the  site  contain  uncorrectable 
spikes  or  clipped  data. 

• One  or  more  components  contain  power  surges  raising  the 
segment  powers  more  than  ten  times  the  power  of  the  preced- 
ing and  following  segments. 

• One  or  more  components  contain  segments  with  zero  power 
(data  drop-outs), 

• One  or  more  components  contain  segment  powers  consistently 
much  higher  than  the  segment  powers  at  other  sites. 

Tne  decision  to  accept  or  reject  a site  is  made  using  the  edit 
printout,  v/hich  lists  segment  powers  for  each  component  of  each  site  edited. 
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Of  the  281  events  in  the  data  base,  unreadable  or  missing  data 
prevented  the  processing  of  29  events.  The  number  of  sites  rejected  from 
beamforming  each  of  the  remaining  252  events  is  summarized  in  Table  IV-2. 
Since  site  1 did  not  become  operational  until  5 May,  there  are  only  232  events 
for  this  site.  Site  3 was  not  operational  during  May. 

We  see  from  this  table  that  site  6 had  the  highest  rejection  rate. 
In  almost  all  cases,  data  from  this  site  was  rejected  due  to  synchronization 
errors  caused  by  the  relaying  of  data  from  site  6 through  site  5 to  the  central 
recording  station.  The  bulk  of  these  synchronization  error  rejections  occurred 
on  the  days  20-23  May  and  30-31  May.  These  errors  reduced  the  site  6 data 
to  meaningless  numbers  which,  when  edited,  caused  many  illegal  gain  codes 
to  be  noted.  When  these  errors  occurred  at  site  6,  the  edit  printout  showed 
much  higher  segment  powers  than  the  other  sites,  with  some  data  drop-outs. 

Sites  1,  2,  and  4 had  approximately  the  same  rejection  rate. 

The  reasons  for  rejecting  these  sites  were  most  commonly  uncorrectable 
spikes  and  power  surges.  Sites  5 and  7 had  the  lowest  rejection  rates.  The 
reason  for  rejecting  these  sites  was  most  commonly  that  one  component  con- 
tained segment  powers  consistently  much  higher  than  the  segment  powers  at 
other  sites. 

C.  SUMMARY  OF  OBSERVATIONS  ON  DATA  QUALITY 

The  remarks  made  to  this  point  on  the  quality  of  the  ILPA  data 
may  lead  one  to  believe  that  the  data  quality  is  poor.  However,  as  was  men- 
tioned earlier,  we  have  found  ways  of  coping  with  most  of  the  problems  we 
found.  The  data-ahift  error  was  easily  corrected  by  aligning  the  data  on  the 
'IL'  characters.  Except  when  timing  errors  occur  at  the  desired  edit  start 
time,  these  errors  do  not  interfere  with  the  running  of  the  edit  program. 

When  timing  errors  do  occur  at  the  edit  start  time,  we  change  the  edit  start 
time  to  a time  a few  minutes  later.  Since  we  always  edit  a noise  gate  before 
the  desired  signal,  this  change  in  the  edit  start  time  does  not  affect  the  signal. 
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The  synchronization  errors  affecting  site  6 are  handled  by 
rejecting  this  site  when  beamforming.  The  only  effect  this  has  is  to  decrease 
the  signal -to-noise  ratio  gain  which  the  beamforming  provides. 

Overall,  we  judge  the  ILPA  data  quality  to  be  fairly  good.  Of 
the  281  events  of  our  data  base,  84.  0 percent  were  successfully  processed, 

1.  8 percent  were  lost  due  to  no  data  being  recorded,  5.  3 percent  were  lost  to 
uncorrectable  malfunctions,  and  8.  9 percent  were  lost  to  unreadable  data. 

The  uncorrectable  malfunctions  were  predominantly  data  spikes  caused  by  bad 
data  associated  with  timing  errors. 


SECTION  V 
SIGNAL  ANALYSIS 


A.  SIGNAL -TO -NOISE  RATIO  GAINS  DUE  TO  BEAMFORMING 

When  the  data  recorded  at  the  individual  sites  of  an  array  are 
formed  into  beams,  the  signal-to-noise  ratio  of  each  component  is  increased 
due  to  suppression  of  noise,  In  the  ideal  case,  the  noise  is  purely  random  and 
is  suppressed  by  a factor  approximately  equal  to  the  square  root  of  the  number 
of  sites  used  in  the  beamforming  process.  In  practice,  the  noise  is  composed 
of  a random  component  and  a propagating  non- random  component.  This  propa- 
gating component  is  suppressed  to  a lesser  degree  than  the  random  component, 
the  amount  of  suppression  depending  on  how  far  off  the  beamforming  azimuth 
its  azimuth  lies.  Also,  the  beamforming  process  suppresses  the  signal  to 
some  extent.  This  is  dependent  on  how  accurate  the  computed  time  delays 
used  to  time-align  the  individual  traces  are.  In  particular,  at  some  point 
close  to  the  array,  the  plane -wave  assumption  used  to  compute  these  time  de- 
lays must  break  down. 

To  obtain  a first  estimate  of  the  signal-to-noise  ratio  gains  we 
may  expect  from  the  beamforming  process,  we  selected  a small  suite  of  events 
which  were  detected  on  both  the  reference  site  and  beam  traces  and  contained 
only  noise  in  the  time  gate  immediately  preceding  the  signal  arrival  time.  The 
signal-to-noise  ratio  for  each  component  of  the  reference  site  and  beam  traces 
were  then  computed  using  the  equation: 

S/N  (SB)  = 20.  * LOG,„  -er°~P“k  a">Pllt"d<’ 

' ' 10  RMS  noise 
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where  "7,fero-to-peak  amplitude"  is  the  amplitude  of  the  largest  peak  of  the 
signal  waveform  and  "RMS  noise"  is  measured  in  the  time  gate  immediately 
preceding  the  signal  arrival.  The  gain  due  to  beamforming  is  then  simply  the 
difference  between  the  beam  signal-to -noise  ratio  and  the  reference  site  signai- 
to-noise  ratio. 

The  results  in  Table  V-l  are  grouped  on  the  basis  of  epicentral 
distance.  The  values  in  the  column  headed  "optimum  gain"  were  computed 
from  the  following  equation: 

Optimum  Gain  {dB)  = 20.  >!<  LOG^  ^number  of  sites 

where  "number  of  sites"  is  the  average  number  of  sites  used  in  beamforming. 
We  see  in  this  table  that  beamforming  gains  for  the  vertical  and  radial  com- 
ponents of  events  with  epicentral  distances  less  than  ten  degrees  are  very  low. 
This  indicates  that  the  plane-wave  assumption  used  in  beamforming  to  compute 
time  delays  fails  for  events  with  epicenters  less  than  ten  degrees  from  the 
array.  The  high  gain  value  for  the  corresponding  transverse  component  is 
not  yet  understood. 

The  mean  gains  for  the  other  ranges  of  epicentral  distances 
remain  fairly  constant.  This  implies  that  the  plane-wave  assumption  holds 
for  events  with  epicentral  distances  greater  than  ten  degrees. 

Comparing  the  mean  gains  in  Table  V-l  with  the  corresponding 
optimum  gains,  we  see  that  in  general  the  mean  gains  are  lower  than  the  opti- 
mum gains.  This  implies  that  some  of  the  noise  1b  propagating,  since,  as 
was  described  earlier,  propagating  noise  is  suppressed  by  beamforming  to  a 
lesser  degree  than  is  random  noise. 

An  interesting  feature  of  the  data  in  Table  V-l  is  that  the  rad- 
ial component  gains  are  lower  than  the  vertical  component  gains.  In  Table 
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TABLE  V-l 

SIGNAL-TO- NOISE  RATIO  GAINS  DUE  TO  BEAMFORMING 


Number  of 
Samples 

Average 
Number 
Of  Sites 

dB  Gains 

Optimum 

Gain 

A Range 

V 

T 

R 

o 

o 

1 

o 

o 

4 

4.  50 

-0.29 

6.  05 

-2.  93 

6.53 

I 0°-20° 

13 

5.69 

7.  67 

4.  50 

3.85 

7.55 

o 

o 

1 

o 

o 

fM 

10 

5.  30 

6.40 

2.  93 

1.  50 

7.24 

>40° 

12 

5.42 

4.  33 

2.76 

3.  95 

7.  34 

t—t 

0 
o 

1 

03 

O 

O 

35 

5.49 

6.  17 

3.70 

3.46 

7.  39 
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V-2  we  have  separated  the  signal-to-noise  ratio  gains  into  the  signal-to-noise 
gain  due  to  RMS  noise  suppression  and  the  signal-to-noise  ratio  loss  due  to 
peak  signal  suppression.  From  the  data  in  this  table  we  see  that  the  differ- 
ence in  signal-to-noise  ratio  gain  between  the  vertical  and  radial  components 
is  due  to  both  lower  RMS  noise  suppression  and  greater  signal  suppression  on 
the  radial  component,  with  the  greater  signal  suppression  being  the  dominant 
factor.  The  lower  RMS  noise  suppression  on  the  radial  component  in  con- 
junction with  lower  RMS  noise  suppression  on  the  transverse  component  rela- 
tive to  the  vertical  component  implies  that  there  is  more  propagating  noise  on 
the  horizontal  components  than  on  the  vertical  components.  We  do  not  as  yet 
understand  the  greater  peak  signal  suppression  on  the  radial  components. 

B.  SHORT-PERIOD  ILPA  DETECTION  CAPABILITY 

During  the  current  contract  period  there  was  no  plan  to  process 
ILPA  short-period  data.  However,  we  were  able  to  review  develocorder 
films  of  ILPA  short-period  data  for  the  first  ten  days  of  May  1976.  This 
period  covered  events  1 to  45  of  our  data  base.  The  detection  results  pre- 
sented here  are  based  on  the  short-period  vertical  component  of  site  7,  (This 
is  the  short-period  component  which  is  recorded  on  the  ILPA  satellite  format 
tapes.  Thus,  our  results  are  indicative  of  the  detection  capability  of  the  un- 
filtered data  available  to  us  on  these  tapes. ) The  detection  criteria  for  these 
data  are: 

• Presence  of  an  impulsive  waveform  occurring  within  + 20  sec- 
onds of  the  predicted  P-wave  arrival  time. 

• Waveform  under  consideration  is  at  least  12  dB  above  the  pre- 
ceding noise. 

The  results  of  this  very  brief  analysis  are  presented  in  Figure 
V-l.  The  upper  portion  of  this  figure  shows  the  distribution  of  detected  and 
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TABLE  V-2 

RMS  NOISE  SUPPRESSION  AND  PEAK  SIGNAL  SUPPRESSION 
DUE  TO  BEAMFORMING 


RMS  Noise  Suppression 
Due  To  Beamforming  (dB) 

Peak  Signal  Suppression 
Due  To  Beamforming  (dB) 

Resultant  Change  In 
Signal-To-Noise  Ratio  (dB) 


A Range 

V 

0 
o 

1 

o 

o 

5.  56 

0 
o 

1 

CD 

o 

o 

6.79 

0 
o 

1 

o 

o 

5.  85 

o 

o 

ao 

1 

o 

o 

0.  62 

o 

o 

1 

c 

o 

-0.  29 

0 
o 

1 

00 

o 

o 
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non-detected  events  as  a function  of  bodywave  magnitude  (m^).  The  l°wer 
portion  of  this  figure  shows  the  maximum  likelihood  curve  fitted  to  these 
detection  statistics  (Ringdal,  1974). 

The  subset  of  events  used  in  this  computation  of  detection  capa- 
bility has  a mean  epicentral  distance  of  37.  1 degrees  with  a standard  deviation 
of  23.  8 degrees.  These  values  are  essentially  the  same  aB  those  for  entire 
data  base.  From  Figure  V-l,  we  see  that  these  events  give  us  an  ILPA  short- 
period  BO  percent  detection  threshold  estimate  of  4,08  m^  units. 

C.  INDIRECT  ESTIMATES  OF  LONG-PERIOD  DETECTION  CAPABILITY 

It  is  possible  to  derive  detection  threshold  magnitudes  from 
ambient  noise  levels,  since  detection  of  a seismic  event  depends  on  the  signal  - 
to-noise  ratio  at  the  recording  station.  Unger  (Unger,  1974)  developed  the 
theoretical  background  for  this  method  and  tested  it  on  Very  Long  Period  Ex- 
periment data.  In  this  method,  it  ie  assumed  that  an  event  can  be  detected 
when  its  maximum  amplitude  exceeds  that  of  the  surrounding  noise  by  a cer- 
tain margin.  The  detection  capability  estimation  algorithm  which  Unger 
develops  is: 


M _ = MEAN  LOG  AMP  - log,.  T *G(T  ) + logind  + d(T  ) - b + C 

b 50  °10  o o °10  o o 

where 


Ms  50 

MEAN  LOG  AMP  = 


G<T0) 


the  50  percent  surface  wave  detection  threshold, 
the  mean  of  the  logarithms  of  the  maximum  peak- 
to-peak  seismometer  output  noise  amplitudes, 
the  geometric  mean  of  the  period  of  the  maximum 
signal  amplitude, 

the  instrument  response  correction  for  period  T , 
taken  from  instrument  response  curves  supplied 
to  us  by  Wayne  Ellis  of  Texas  Instruments  Incor- 
porated in  Dallas,  Texas, 
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Aq  = the  geometric  mean  of  the  signal  epicentral  dis- 

tances, 

d(T  ) = the  station  magnitude  difference  due  to  period  T , 

b = the  mean  station  bias, 

C =1.12+  the  detection  criterion  margin. 

We  used  17  noise  samples  from  May  to  compute  MEAN  LOG 
AMP.  The  parameter  was  selected  to  be  20  seconds,  since  this  period 
was  more  often  observed  in  detected  signals  than  either  30  seconds  or  40  sec- 
onds. G(Tq)  was  then  the  instrument  response  correction  at  20  seconds  picked 

from  the  long-period  instrument  response  curve.  The  parameter  A was 

o 

varied  from  10  to  100  in  10  increments  to  give  us  a table  of  M^  ^ versus 
epicentral  distance.  The  parameter  d(TQ)  was  picked  from  the  plot  of  mag- 
nitude difference  versus  period  found  in  Unger's  report  (Unger,  1974).  For  a 
period  of  20  seconds  and  a continental  path,  this  parameter  has  the  value  of 
-0.01.  Since  we  had  no  previous  knowledge  of  the  mean  station  bias,  b was 
assumed  to  be  zero,  The  detection  criterion  was  set  at  a factor  of  two,  i.  e.  , 
the  maximum  signal  amplitudes  must  be  twice  the  maximum  noise  amplitudes, 
giving  us  C=  1.12+  0.  301  = 1.421.  The  results  of  these  computations  are 
given  in  Table  V-3. 

D.  DIRECT  ESTIMATES  OF  LONG-PERIOD  DETECTION  CAPABILITY 

In  this  subsection,  we  shall  examine  the  direct  estimates  of 
detection  capability  of  the  reference  site  and  beam  data.  The  criteria  which 
determine  whether  an  event  was  detected  are; 

• The  presence  of  signal  dispersion  in  the  signal  gate, 

• A peak  in  the  dispersed  wavetrain  3 dB  or  more  above  any 
peak  outside  the  dispersed  wavetrain  and  inside  a 20-minute 
time  gate  centered  at  the  expected  peak  occurrence  time  , 
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TABLE  V-3 


INDIRECT  ESTIMATES  OF  LONG -PERIOD  DETECTION  CAPABILITY 

VERTICAL  COMPONENT 


Reference  Site 

Beam 

Ao 

U 

M 50 

s 

m,  50 
b 

M 50 
s 

m,  50 
b 

10° 

2.  78 

3.  91 

2.  50 

3.74 

20° 

3.  08 

4.  08 

2.  80 

3.92 

30° 

3.  26 

4.  19 

2.  98 

4.  02 

o 

o 

3.  38 

4.26 

3.  10 

4.09 

50° 

3.  48 

4.  31 

3.  20 

4.  15 

60° 

3.  56 

4.  36 

3.28 

4.20 

o 

o 

3.  63 

4.  40 

3.  35 

4.24 

oo 

o 

o 

3.  68 

4.43 

3.  40 

4.27 

o 

o 

CT- 

3.74 

4.47 

3.  46 

4.  30 

100° 

3.78 

4.49 

3.  50 

4.  3 3 

m,  50  is  calculated  from  M s 1,72  m,  - 3.  94 
b p b 


(Table  V-5  vertical  component  20-second  period) 
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• The  onset  of  the  signal  occurs  within  +.180  seconds  of  the  ex- 

pected signal  onset  time. 

Of  the  2.52  events  which  were  successfully  processed,  we  found 
no  events  which  we  designated  as  detected  which  did  not  fulfill  the  first  two 
criteria.  Occasionally,  the  third  criterion  was  not  fulfilled,  the  observed 
waveform  arriving  later  than  expected.  In  such  cases,  the  event  lists  were 
checked  to  ensure  that  the  detected  waveform  was  not  due  to  some  other 
event.  If  no  other  event  could  be  found  whose  surface  waves  would  arrive  at 
that  time,  the  event  under  consideration  was  called  a detection. 

We  first  calculated  the  detection  capability  of  the  array  using 
all  events  of  the  data  baoe.  In  this  case,  mixed  events,  events  containing 
malfunctions,  events  for  which  no  data  was  recorded,  and  events  for  which 
the  data  was  unreadable  were  ail  counted  as  non-detections.  The  results  for 
this  case  are  shown  in  Figures  V-2  to  V-4,  The  upper  portion  of  each  figure 
shows  the  distribution  of  detected  and  non-detected  events  as  a function  of 
bodywave  magnitude  (m^).  The  Lower  portion  of  each  figure  shows  the  maxi- 
mum likelihood  curve  fitted  to  these  detection  statistics.  In  the  figure, 
"MB50"  denotes  the  BO  percent  detection  threshold,  "MB90"  denotes  the  90 
percent  detection  threshold,  "SIGMA”  is  related  to  the  slope  of  the  maximum 
likelihood  curve,  and  "RHO"  denotes  the  quality  of  the  results  (Ringdal,  1974). 

These  figures  tell  us  the  absolute  probability  of  detecting  an 
event  with  given  from  Eurasia.  For  example,  the  50  percent  detection 

probability  is  at  m^  equal  to  4.  55  for  Rayleigh  waves  and  4.48  for  Love  waves. 

We  next  consider  the  case  where  only  events  for  which  a de- 
tection/non-detection decision  could  be  made  were  included  in  the  detection 
statistics.  This  gives  us  the  conditional  detection  probability  curves  shown 
in  Figure  V-5  for  reference  site  data  and  in  Figure  V-6  for  beam  data.  (The 
term  'conditional'  means  that  the  probability  curves  are  computed  under  the 
condition  that  a clear  detection/non-detection  decision  he  possible  for  each 
event  considered.  ) We  do  not  show  the  corresponding  plots  for  the  tranBverBe 
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ILPA  BEAM  VERTICAL  COMPONENT  DETECTION  STATISTICS 
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and  radial  components  because  the  results  are  essentially  the  same.  The  de- 
tection capability  is  summarized  in  Table  V-4.  We  note  that  the  indirect  esti- 
mates were  taken  from  Table  V-3  using  the  mean  epicentral  distance  for  the 
data  base  (34,5°).  These  indirect  estimates  are  conditional  detection  capa- 
bility estimates  since  the  assumption  that  there  are  no  malfunctions,  mixed 
signals,  missing  data,  or  unreadable  data  is  implicit  in  their  computation. 

The  following  points  concerning  the  data  in  Table  V-4  can  be 

made: 

• Malfunctions,  mixed  signals,  miBsing  data,  and  unreadable 
data  cost  us  0.45  m^  units  in  detection  capability. 

• The  indirect  detection  estimates  are  in  close  agreement  with 
the  direct  estimates. 

• Beamforming  lowers  the  50  percent  detection  threshold  by  ap- 
proximately 0.  15  units  on  all  components, 

Although  the  data  base  was  too  small  to  regionalize  the  detec- 
tion capability,  it  was  possible  to  approximate  this  by  separating  the  events 
into  three  groupB  based  on  epicentral  distance.  These  groups  are  near  field 
(0°2JA<20°),  near  teleseismic  (20°5: A<  50°),  and  far  teleseismic  (50°S»A  <80°). 
The  maximum  likelihood  detection  curves  for  these  three  cases  are  Bhown  in 
Figures  V-7  to  V-9.  Again,  these  curves  give  conditional  detection  probabil- 
ities, since  all  mixed  events,  events  with  malfunctions,  events  with  missing 
data,  and  events  with  unreadable  data  have  been  removed  from  the  detection 
statistics.  As  before,  we  Bhow  only  the  vertical  component  results,  since  the 
results  for  the  horizontal  components  were  essentially  the  same. 

Up  to  this  point  we  have  diBcussed  detection  capability  in  terms 
of  individual  components.  We  will  now  lpok  at  the  combined  detection  capabil- 
ity of  the  components.  (Again,  we  will  use  conditional  detection  probabilities, 
i.  e.  , all  mixed  events,  events  for  which  no  data  was  recorded,  events  which 
contained  malfunctions,  and  events  for  which  the  data  was  unreadable  have 
been  removed  from  the  detection  statistics.  ) In  Figure  V-10  we  show  the 
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TABLiE  V -4 

SUMMARY  OF  DETECTION  CAPABILITY  INFORMATION 
50  PERCENT  BODYWAVE  MAGNITUDE  (mb)  DETECTION  CAPABILITY 


Component 

Absolute 

Detection 

Capability 

Conditional  Detection  Capability 

Indirect  Estimates 

Direct  Estimates 

Reference 

Site 

Beam 

Reference 

Site 

Beam 

Beam 

Vertical 

4.  55 

4.  22 

4.  05 

4.  25 

4.  09 

Transverse 

4.  48 

* 

♦ 

4.  16 

4.  04 

Radial 

4.  55 

* 

* 

4.  24 

4,  09 

Not  Measured 
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IUPA  beam  vertical  component  detection  statistics 

0°  £ A < 20° 
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FIGURE  V-10 

ILPA  BEAM  DETECTION  STATISTICS:  DETECTION  ON  ONE 
COMPONENT  REQUIRED  TO  DECLARE  AN  EVENT  DETECTION 
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beam  detection  statistics  for  the  case  where  an  event  is  declared  to  be  detect- 
ed if  it  is  detected  on  one  or  more  component.  The  50  percent  detection  thresh- 
old for  this  case  is  at  m^=  4.  01,  We  might  call  this  case  our  optimistic  detec- 
tion capability  estimate. 

In  Figure  V-ll  we  Bhow  the  beam  detection  statistic s for  the 
case  where  an  event  is  declared  to  be  detected  only  if  it  is  detected  on  all  three 
components.  The  50  percent  detection  threshold  for  this  case  is  at  m^=  4.  13. 
We  might  call  this  case  our  conservative  detection  capability  estimate. 

E.  M - m,  RELATIONSHIP  AT  ILPA 

We  cannot  estimate  the  discrimination  capability  of  ILPA  since 
our  data  base  does  not  contain  any  presumed  explosions.  We  can,  however, 
make  a start  on  this  portion  of  the  ILPA  evaluation  by  describing  the  earth- 
quake M - m.  relationships  at  this  array, 
s u 

Surface  wave  magnitudes  (M  ) were  computed  using  the  equa- 

s 

tions 

Ms  = L°g10  [ S.  F.  "*“G  * Q * T ] + L°g10A  + 1,12 

where 

A = Peak-to-peak  amplitude  measured  in  inches  on  the  plot  of  the 
event, 

S.  F.  = The  plot  scale  factor  in  inches  per  computer  count, 

G = The  instrument  response  correction  factor, 

Q - The  quantization  rate  (20.951  cc/m //), 

T = The  period  at  which  the  amplitude  A was  measured, 

A = The  epicentral  distance  of  the  event. 
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ILPA  BEAM  DETECTION  STATISTICS:  DETECTION  ON  ALL 
COMPONENTS  REQUIRED  TO  DECLARE  AN  EVENT  DETECTION 
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The  normalized  instrument  response  curve  from  which  G was 
taken  is  shown  in  Figure  V-12.  This  curve  was  created  by  averaging  the  indi- 
vidual site  instrument  response  curves  supplied  by  Wayne  Ellis  of  Texas  In- 
struments Incorporated  in  Dallas.  (One  normalized  instrument  response 
curve  for  each  site  was  supplied.  Since  this  was  the  only  information  we  had, 
we  used  this  curve  for  all  three  components.  ) The  vertical  bars  on  the  curve 
show  plus  or  minus  one  standard  deviation  at  each  of  the  measurement  points. 


Wherever  possible,  was  measured  at  periods  of  20,  30,  and 
40  seconds.  Figures  V-13  to  V-15  show  the  M -m,  data  with  M measured  at 
20  seconds  on  the  vertical,  transverse,  and  radial  beams,  respectively.  The 
fit  to  the  data  points  shown  in  each  figure  treats  neither  variable  as  depend- 
ent, minimizing  the  perpendicular  distance  from  the  line  to  the  data  points. 


In  Table  V-5  we  list  the  slope  and  intercept  of  these  M^-m^ 
fits,  where  a and  b are  the  coefficients  in  the  equation: 


Mg  = a m^  + b. 


We  used  this  information  to  construct  Figure  V-16,  showing  the  fits 

at  20,  30,  and  40  seconds  for  ec,ch  component  of  motion.  These  plots  show 

that,  on  the  average,  the  30  second  M is  lower  than  the  20  second  M and 

0 0 

the  40  second  M is  lower  than  both  the  20  second  and  30  second  M . 

s s 

At  this  point  we  are  in  need  for  depth  information  for  the  earth- 
quakes and  a suite  of  presumed  explosions  to  carry  out  the  evaluation  of  the 
ILPA  discrimination  capability.  The  depth  information  will  be  used  to  elimi- 
nate deep  events  from  the  data  base.  These  events  should  be  removed,  since 
they  tend  to  generate  lower  surface  waves  than  do  the  shallow  events.  Thus, 
their  removal  will  lower  the  variance  of  the  earthquake  population  about  the 
Mg-m^  earthquake  fit.  The  presumed  explosion  data  would  then  be 

plotted,  and  we  could  finally  determine  the  quality  of  this  discriminant. 
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TABLE  V-5 

mb  FIT  SLOPE  AND  INTERCEPT  VALUES 


Period 

(sec) 

Component 

a 

b 

tf2 

Center  of  Mass 

n 

mb 

M 

8 

Vertical 

1.72 

-3.  94 

0.  10 

4.  35 

3.  55 

62 

20 

Transverse 

1.68 

-3.61 

0.  11 

4.  33 

3.  66 

66 

Radial 

1.64 

-3.59 

0.  09 

4.  32 

3.  50 

64 

Vertical 

1.72 

-4.  15 

0.  10 

4.  34 

3.  30 

61 

30 

Transverse 

1. 59 

-3.42 

b.  13 

4.  33 

3.  47 

61 

Radial 

1.71 

-4.  18 

0.  10 

4.  38 

3.  29 

61 

Vertical 

1.51 

-3.  73 

0.  10 

4.46 

2.  99 

29 

40 

Transverse 

1.51 

-3.  57 

0.  09 

4.  40 

3.  09 

36 

Radial 

1.79 

-4.  98 

0.  11 

4.  42 

2.  92 

24 

where  M = am.  + b 
s b 
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SECTION  VI 
CONCLUSIONS 


In  thi8  section,  we  will  summarize  the  results  determined  in 
this  first  evaluation  of  the  Iranian  Long  Period  Array.  The  major  points  are: 

• The  ILPA  data  quality  is  judged  t:o  be  fairly  good  despite  the 
various  problems  encountered  in  reading  the  data  tapes,  Of 
the  281  events  in  the  data  base,  84.0  percent  were  success- 
fully processed,  1.  8 percent  were  lost  due  to  no  data  being 
recorded,  5.  3 percent  were  loot  to  uucorrectable  malfunc- 
tions, and  8,  9 percent  were  lost  to  unreadable  data. 

• The  highest  gain  in  signal -to -noise  ratio  due  to  boauiforming 
was  6.  17  dB  on  the  vertical  component.  Gains  on  the  hori- 
zontal components  averaged  between  3 and  4 dB. 

a The  ILPA  short-period  50  percent  detection  threshold  was 

estimated  to  be  4.08  cn^  units.  Detection  statistics  used  in 
making  this  estimation  were  obtained  by  visually  reviewing 
develocorder  films  for  events  1 to  45. 

• The  absolute  50  percent  detection  threshold  for  ILPA  beam 
data  is  at  m^  = 4.  55  for  Rayleigh  waves  and  = 4.  48  for 
Love  waves.  Absolute  detection  thresholds  were  arrived  at 
by  counting  all  mixed  events,  events  for  which  no  data  was 
available,  events  containing  malfunctions,  and  everts  for 
which  the  data  was  unreadable  as  non-detections. 
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• The  conditional  50  percent  detection  threshold  for  ILiPA  beam 
data  is  at  4.  09  for  Rayleigh  waves  and  4,  04  for  Love 
waves.  Conditional  detection  thresholds  were  arrived  at  by 
including  in  the  detection  statistics  only  those  events  for  which 
a detection/ non-detection  decision  could  be  made. 

• Beamforming  lowered  the  conditional  50  percent  detection 
threshold  by  approximately  0.  15  units. 

• Indirect  estimates  of  the  50  percent  detection  threshold  made 
from  noise  samples  agreed  quite  closely  with  the  above  de- 
scribed direct  estimates. 

• Combining  the  beam  detection  statistics  of  the  individual  com- 
ponents resulted  in  a conditional  50  percent  detection  thresh- 
old of  m.  = 4.  01  for  the  case  where  an  event  was  declared  de- 

b 

tected  if  it  was  detected  on  one  or  more  components.  This 
conditional  50  percent  detection  threshold  is  4.  13  for  the 
case  where  the  event  was  declared  detected  only  if  it  was  de- 
tected on  all  three  components. 

• m^  fits  were  computed  at  periods  of  20,  30,  and  40  sec- 
onds for  earthquake  data  only,  since  the  data  base  did  not  con- 
tain any  presumed  explosions.  For  this  data  base,  the  surface 
wave  magnitude  decreased  with  increasing  period. 

The  major  areas  which  must  be  investigated  in  the  future  to 
complete  the  evaluation  of  the  Iranian  Long  Period  Array  are  as  follows: 

• Noise  analysis  - A daily  sampling  cf  the  noise  field  is  need- 
ed to  provide  us  with  estimates  of  RMS  noise  levels,  spectral 
content  cf  the  noise,  and  noise  coherency. 
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Signal  analysis  - The  data  base  must  be  greatly  expanded  so 
that  regional  detection  capability  can  be  estimated.  Also,  the 
work  on  signal-to-noise  ratio  gains  due  to  beamforming  should 
be  continued  and  expanded. 

Discrimination  capability  - Depth  Information  for  the  earth- 
quakes and  a suite  of  presumed  explosions  must  be  obtained  to 
determine  the  ILPA  discrimination  capability. 
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APPENDIX  A 
THE  DATA  BASE 


In  this  appendix  we  present  the  parameters  describing  each  of 
the  events  of  the  data  base.  The  column  headed  "EVNO"  gives  the  unique 
number  assigned  to  each  event.  The  column  headed  "DATE"  gives  the  month, 
day,  and  year  of  occurrence  of  the  event.  The  column  headed  "TIME"  giveB 
the  event  origin  time.  The  columns  headed  "LAT"  and  "LONG"  give  the  lati- 
tude and  longitude  of  the  event  epicenter,  where  a positive  value  indicates 
north  latitude  or  east  longitude  (as  appropriate)  and  a negative  value  indicates 
south  latitude  or  west  longitude.  The  column  headed  "MB"  gives  the  NORSAR 
value  of  the  event  bodywave  magnitude.  The  column  headed  "LOCATION" 
gives  the  general  seismic  region  in  which  the  event  occurred.  Finally,  the 
column  headed  "QUALITY"  gives  the  NORSAR  quality  rating  of  the  event 
parameters,  where 

1 = good  to  excellent 

2 a fair  to  good 

3 = poor  to  fair. 

We  note  that  all  the  event  parameters  listed  wer.-  . ■ en  from  the  NORSAR 
event  lists. 
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